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Abstract 

The rheology of granular materials near an interface is investigated through proton magnetic resonance imaging. A 
new cylinder shear apparatus has been inserted in the MRI device, which allows the control of the radial confining 
pressure exerted by the outer wall on the grains and the measurement of the torque on the inner shearing cylinder. 
A multi-layer velocimetry sequence has been developed for the simultaneous measurement of velocity profiles in 
different sample zones, while the measurement of the solid fraction profile is based on static imaging of the sample. 
This study describes the influence of the roughness of the shearing interface and of the transverse confining walls on 
the granular interface rheology. 
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1. Introduction 

The interaction of a granular material with a solid in- 
terface is of interest in various engineering problems, 
such as industrial conducts geotechnics ylEL and 
in geophysical situations, such as tectonophysics J3l or 
gravity flows Eh- 

At the immediate vicinity of the shearing interface, 
a thin granular layer, where the shear and dilation is lo- 
calized, plays a significant role in the stress transmission 
between the solid interface and the bulk granular mate- 
rial. This rheology is influenced by the roughness of the 
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In this paper, we focus our attention on the annular 
(Couette) shear geometry, where the material is con- 
fined between two cylinders and sheared by the ro- 
tation of the inner rough one (see Jl3ll for a recent 
review). This geometry has been used to measure 
the rheological properties of granular materials, both 
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ever, the visualization of the granular interface is usu- 
ally limited to the upper (free surface) or bottom layers 
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der vibration J40l |4lj, |42j SlIIl SH we have 
used MRI to measure the granular rheology (velocity 
and solid fraction profiles) well inside the sample. 

Sec.|2]is devoted to the description of a new annular 
shear cell, specially designed to be inserted in a MRI 
device. Based on a geotechnical cylinder shear appara- 
tus (Hi! |H 0, El EH, its originality relies on 
the control of the radial confining pressure exerted by 
the outer wall on the grains and on the measurement of 
the torque on the inner shearing cylinder. 

Sec. [3] explains the multi-layer MRI velocimetry. 
MRI velocity measurements were performed using 
a spin-warp / phase encoding technique previously 
adapted from 115 211 and used for routine liquid rheol ogy 
in annular Couette cells (see for instance lf53"ll54 , 55, 5a 



5710 . It was here further modified on purpose of quasi- 



simultaneous assessment of different regions in the sam- 
ple. 

Sec.|4]describes the measurement of the velocity and 
solid fraction profiles, as well as the shear stress on the 
shearing wall, from which we deduce the influence of 
the roughness of the shearing interface and of the trans- 
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verse confining walls on the granular interface rheology. 
Our study is restricted to the quasi-static regime (small 
shear velocity and/or high confining pressure). 

Preliminary and complementary results are presented 
in 0|. 

2. Annular shear cell 

The annular shear cell is inserted in the MRI de- 
vice through an external support. The figure Q] shows a 
schematic view of the cell, which components are made 
ofPMMA. 
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Figure 1: Schematic view of the shear cell. 



The sample has a hollow cylinder shape. The gran- 
ular material is confined between two fixed horizontal 
(bottom and upper) plates (height H — 10 cm), an in- 
ternal rotating cylinder (radius Rj = 3 cm) and an ex- 
ternal latex homemade dip molded membrane (radius 
R e — 6 cm) filled with water. A pressure-volume con- 
troller (GDS®) applies a radial confining pressure P to 
the sample, in the range 0-15 kPa. An optical fiber 
sensor (FOP MEMS 1000 kPa - Fiso Technologies®) 
precisely measures this pressure close to the membrane. 

The internal rotating cylinder is guided by two jour- 
nal bearings set in the external support, while the main 
cell (membrane, bottom and upper plates) is only con- 
nected to the external support through the torque sensor, 
which measures the whole torsion effort. 

An aluminum-alloy torque sensor was specially de- 
signed, based on resistive strain gages, in the range 



+ 10 Nm. This measurement was not possible during 
MRI experiment, but when displacing the cell bellow 
the MRI device. This prototype is a first step toward the 
realization of a torque sensor working inside the MRI 
device. 

Depending on the way upper and bottom plates are 
mounted, it is possible to measure the whole torsion ef- 
fort or only the fraction transmitted to the lateral mem- 
brane (the difference between those two measurements 
provides the torque transmitted by the horizontal walls). 

The cell is connected to the transmission axis of a 
rheometer previo usly designed to be inserted within the 
MRI facility JH ES El H, through a gearbox. This 
two-stage, timing-belt and pulley system is placed close 
to the cell and far from the motor. Its reduction factor of 
10 provides a rotation range 1 /600 < £1 < 1/6 RPS. 

This configuration allows to place down the cell (out 
of the MRI tunnel) during the sample preparation and 
torque measurement and then move it up to the observa- 
tion position. 

The complete cell has a total diameter of * 19.5 cm 
and a total height, without gearbox, of 28 cm, which fits 
inside the RF coil. 



3. Multi-layer MRI velocimetry 

MRI experiments were carried out on a Bruker 
Biospec 24/80 MRI facility operating at 0.5 T (21 MHz 
proton frequency). The MRI magnet is a vertical su- 
perconducting prototype (Magnex Scientific, Oxford), 
with a 40 cm bore. These characteristics are particu- 
larly suited for the study of large and inhomogeneous 
samples, exhibiting strong internal susceptibility con- 
trasts. The magnet is equipped with a birdcage RF coil 
(height: 20 cm, inner diameter: 20 cm, hard n/2 pulse 
duration: 100jus), and a BGA26 shielded gradient sys- 
tem (Bruker), delivering a 0.05 T/m gradient strength 
with a rising time of 500//.?. 

MRI methodolo gy f or radial velocimetry inside the 
cell was that of ll52ll as further modified by lf59fl 
and IrSOll . It is based on a two-pulse spin echo se- 
quence (Fig. [2^), in which the two pulses are made 
space-selective in z and y direction respectively so as to 
virtually cut a beam along one cell diameter (Fig. |2j)). 
A read-out gradient in x directions permits to get af- 
ter Fourier transform of the signal a ID magnetization 
profile along the beam. An additional pair of gradient 
pulses in y direction (in black) induces on the magneti- 
zation profile an x dependent phase shift reading: 



(p(x) = yG6Av y (x), 



(1) 



where y is the gyromagnetic ratio of proton, G the 
strength of these 'velocity' pulses, 5 and A are timing 
components of the sequence as described in Fig. [2] and 
v y (x) is the y velocity component along the selected di- 
ameter. In order to get a velocity profile, one performs 
two MRI measurements with respectively positive and 
negative velocity gradients, and then compares in each 
pixel the phase of the two magnetization profiles. When 
the thickness of the beam in y direction is small enough 
as compared with the cell diameter, such measurement 
may be regarded as a direct measurement of the orthora- 
dial velocity component vg versus the radial coordinate 
r (Fig. Eh). 
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Figure 2: (a) MRI sequence used for velocity measurements, and its 
coherence pathway. Velocity encoding gradients are shown in black, 
(b) Cell scheme and virtual beam cut. 

In this sequence, only the phase of NMR signal bears 
relevant information. It is then safe to use recycling 
times TR equal or even shorter than T\ values. Indeed, 
an incomplete spin-lattice relaxation between two suc- 
cessive sequences may only affect the signal to noise ra- 
tio, but is not prone to introduce any bias, provided that 
the magnetization be correctly cleaned from past solic- 
itations before each sequence. Such cleaning was here 
performed with a saturation module at the beginning of 
the TR delay. 



Thanks to saturation, TR may even be chosen so as 
to optimize the signal to noise ratio of the experiment. 
Let's imagine that, because of experimental constraints, 
one has a fixed delay T exp to perform a velocity mea- 
surement. If some TR delay is used, then it will be 
possible to repeat the sequence N = T exp /TR times. 
Granted that the available magnetization at the begin- 
ning of each sequence scales as 1 - exp(—TR/T\), and 
that the improvement of signal to noise ratio when re- 
peating sequences scales as yfN, the signal to noise ratio 
SNR of the whole measurement process scales as: 

SNR(TR) cell- ex p(-^)} y^f- ( 2 ) 
This quantity is maximum for: 

TR opt = 1.267!, (3) 

and fairly stays above 90% of this maximum in the 
range 0.577^ < TR < 2.597" i . Saturation then turns 
out a useful tool for signal optimization. 

For the assessment of the velocity profiles at differ- 
ent height in z direction, sequences were repeated while 
systematically shifting the vertical position of the beam. 
In order to study unsteady granular flows, it was impor- 
tant to make these multi-layer measurements as simul- 
taneous as possible. 

First of all, looking back at the NMR sequence, the 
measurement of one beam at a given position clearly af- 
fects all the magnetization of a larger region in the cell, 
composed of the sum of the two slices selected by each 
RF pulse. This, together with the use of the saturation 
module, prevents considering measurements at different 
z positions as independent: it is not possible to run the 
sequence at one height immediately after running the 
sequence at another height, even if the two associated 
beams are well separated. Real simultaneous measure- 
ments of velocity in multiple layers were then not pos- 
sible. 

We proceeded instead as follows (Fig [3j. The whole 
measurement consisted in two nested loops. In the first 
loop the vertical location of the beam was shifted af- 
ter each single sequence, so as to go through all neces- 
sary beam positions. The procedure was then just re- 
peated with inverted velocity gradient strength. This 
first loop was inserted in a second one consisting in stan- 
dard phase cycling and signal accumulation. One could 
get this way velocity profiles at different height mea- 
sured over exactly the same period of time. 

In the present work, the thickness t of the beam in y 
direction was 10 mm, and the recycling time TR was 
taken as 100 ms. The typical echo-time was about 
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Figure 3: Organization of MRI velocimetry sequences for the quasi- 
simultaneous assessment ofn layers in the cell. 



The granular material is poured inside the cell with 
a funnel from the top and then levelled off (see Figj4j. 
After closing the cell, the sample is confined by the ap- 
plication of the radial pressure through the membrane. 
Then, the material is pre-sheared (at least 15 rotations of 
the cylinder, corresponding to « 2800 mm of tangential 
displacement of the wall) to drive it in a steady-state. 
Then data acquisition starts. 




20 ms. Strong enough velocity gradients were here suf- 
ficient to select the proper coherence pathway. Our MRI 
system did not suffer any visible acquisition quadrature 
defect, so that the phase cycling procedure simply con- 
sisted here in a base-line correction, and could be com- 
pleted in two steps. The minimal value of the height /; 
of the beams was 1 .5 mm. The maximum relevant value 
was the cell size. Thinner layers are well adapted to 
study regions with a strong vertical gradient of veloc- 
ity, like the zone close to the horizontal walls. Because 
of hardware constraints, the maximum number of layers 
was 10. The total experimental time was 13 minutes for 
three layers of 1.5 mm width, 104 minutes for five lay- 
ers of 0.4 mm width and 208 minutes for five layers of 
0.15 mm width. 

4. Experimental results 

Three roughness levels of the rotating cylinder are 
tested: the first one corresponds to a vertically corru- 
gated surface (triangular shape with thickness equal to 
depth of 1 mm), the second one to a glued seeds sur- 
face and the third one to sandpaper (medium grain - 80). 
The horizontal walls (upper and bottom) are originally 
smooth. To modify their roughness we have sticked 
sandpaper and glass beads with the same diameter as 
the seeds (for practical reasons the beads are glued to 
sandpaper and not directly over the plate surface). 

Mustard seeds are used as model granular material 
(mean diameter d — 1 .5 mm, quite monodisperse, mass 
density of 1200 kg/m 3 ). The relaxation times are T\ « 
100 ms and T-> as 40 ms. 




Figure 4: (Color on line) Sample preparation. 



Tab. Q] summarizes the parametric study (wall rough- 
ness, confining pressure P and velocity at the wall Vg = 
IttRiQ). 
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1 pressure (kPa) 2 velocity (mm/s) 

Table 1: List of velocimetry experiments. 

Internal stresses associated to gravity (of the order of 
0,7 kPa) can be considered to be dominated by the ap- 
plied confining pressures. 



The maximum value of the inertial number 11311. used 
to qualify the granular flow regime, can be estimated 
around 3 ■ 10~ 4 , indicating that the granular material is 
in the quasi-static regime. 

We now discuss the radial profiles of velocity and 
solid fraction. They are presented as a function of r—Ri, 
normalized by the diameter d of the grains. 



4.1. Velocity profiles 

We have first measured the profiles of the orthoradial 
velocity vg(r), normalized by Vg, far from the horizontal 
walls, comparing three layers 15 mm thick separated by 
10 mm (Pig. |5j. 
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Figure 5: Scheme of layers in the central zone. 



Fig|6] shows that vg(R/) < Vg, which signifies that 
there is significant sliding at the wall. Fig|6]also shows 
that there is no influence of the roughness of the hori- 
zontal walls in the central zone. An approximately ex- 
ponential decay of the orthoradial velocity profiles vg(r) 
is measured, with a maximum value at the shearing wall, 
consistently with previous studies 118 , 
ESQ. Qualitatively, this shear localization is ex- 
plained by the strong decrease of the shear stress away 
from the inner wall, while the normal stress remains ap- 
proximately constant. As expected in the quasi-static 
regime, no influence of the velocity at the wall Vg nor of 
the confining pressure P has been observed. 



4.1.1. Influence of the roughness of the horizontal walls 
In a 2D system, the shear stress is entirely transmit- 
ted from the inner to the outer wall by the granular ma- 
terial 10, [T3I]. In a 3D system, the granular material in- 
teracts with the horizontal walls, to which part of the 
torque is transmitted. This changes the stress distri- 
bution within the sample, and may consequently affect 
the 3D velocity profile. This influence of frictional lat- 
eral walls has been evidenced in granular surface flows 
between vertical walls 06 ill , especially when they are 
rough Il62ll . We now wonder if such effect is observed 
close to the horizontal walls. 

The analysis of the influence of the horizontal walls 
was based on a comparison of the profiles of 5 lay- 
ers 1.5 mm thick separated by a distance of 1.75 mm 
(Fig. |7J. The gap of 0.45 mm between the horizontal 
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Figure 6: (Color on line) Velocity profiles (glued seeds cylinder), (a) 
Smooth horizontal walls (layers 1, 2 and 3). (b) Influence of the rough- 
ness of the horizontal walls (layer \): (») smooth, (®) sandpaper, (o) 
glued glass beads. In inserts, the region closer to the shearing cylin- 
der in semi-logarithmic scale. 



wall and the first layer ensured that only moving parti- 
cles were considered in the measure. In the rough case, 
a similar gap was taken from the glued glass beads sur- 
face giving a total thickness of 2.5 mm. 

The horizontal roughness made of glass beads is not 
detectable by MRI, which provided a correct evaluation 
of the velocity of the grains in direct contact with this 
surface. The motion of the grains could be clearly dis- 
tinguished from the motion of the cylinder, given the 
strong velocity gradient. 

Fig. |8] shows a comparison between the profiles in 
the central zone and near the horizontal wall. In the 
case of a smooth horizontal wall, there is no significant 
influence. This result has practical applications, vali- 
dating the measure of displacements of g ranular mate- 
rials through transparent glass walls Il9l I7L 12911 . Con- 
versely, in the case of a rough horizontal wall, a signif- 
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Figure 7: Scheme of layers close to the inferior plate. 



icant decrease is observed very close to the horizontal 
wall, particularly in the first layer. This indicates the 
transmission of shear stress between the granular mate- 
rial and the horizontal walls. However, the perturbation 
remains very localized, in contrast with granular surface 
flows J62|]. The fluctuations of the velocity profiles in 
the first layers close to the horizontal wall is explained 
by the small thickness of those layers, which strongly 
decreases the signal to noise ratio. 



4.1.2. Influence of the cylinder roughness 

Fig. [9] shows that the cylinder roughness has a sig- 
nificant influence on the velocity vg, causing offsets of 
the whole velocity profiles. Sliding, defined as the ra- 
tio between the maximum value of the particle velocity 
and the wall velocity Vg 1E9I 131. decreases from sand- 



paper cylinder to the corrugated one, that is to say as the 
roughness increases. 

4.2. Friction coefficient 

The average shear stress at the inner cylinder is 
simply deduced from the measured torque C as S — 
C l(2nR^H). Consistently with previous studies 163, 5C 




a linear relation between the shear stress S and 
the confining pressure P is observed (Fig. [Tot , from 
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Figure 8: (Color on line) Velocity profiles near the bottom horizontal 
wall (semi-logarithmic scale, glued seeds cylinder), (a) Smooth hor- 
izontal wall (layers described in Fig. f^'j. (b) horizontal wall with 
glued glass beads (layers described in Fig. [2>). (I) layer 1, (•) 
layer 2, ( k) layer 3, (J) layer 4, (+) layer 5, (O) central profile. 



which an effective friction coefficient S/P is deduced 
(considering that the normal stress at the cylinder wall 
is approximately equal to P) for each of the three cylin- 
ders: 0.35 for sandpaper, 0.38 for glued seeds and 0.4 
for corrugation. This increase seems consistent with the 
decrease of the particle sliding shown in Fig. [9] The 
roughness of the horizontal walls does not affect this 
value, which is consistent with the very localized influ- 
ence of this roughness. 

4.3. Solid fraction profile 

In order to measure the influence of the shear on the 
radial distribution of grains, 3D spin density pictures 
(20 layers 5 mm thick in z direction, with a space reso- 
lution of 0.2 mm in both x and y directions) were per- 
formed at rest before and after shearing (at least 200 
rotations, making « 38 m of tangential displacement). 
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Figure 9: (Color on line) Influence of the cylinder roughness on the 
velocity profiles: (M) sandpaper, (») glued seeds and (A) corrugated, 
(a) Linear scale and (h) detail of the region near the shearing cylinder 
in semi-logarithmic scale. 



Fig.fTTh and Fig.fTTb show typical density pictures ob- 
tained before and after shearing respectively. 

For each radial coordinate, grayscale pictures were 
then integrated together in the orthoradial direction. 
Granted that the integrated spin density profile in radial 
direction provides a relative measurement of the solid 
fraction, these solid fractions v(r) are normalized by the 
mean value v,„ obtained before shearing, and shown in 

Fig.[n] 

Before shearing, after the granular material has been 
poured in the cell, the solid fraction is relatively homo- 
geneous, with a certain structuration (v(r) oscillations) 
near the inner wall. After shearing, a much clearer struc- 
turation (already visible on the static images (Fig. [TTJ) 
together with a significant dilation is observed close to 
the shearing cylinder (r - R, ;$ Id). However, the zone 
Id iS r - Ri iS 12c/ does not reveal significant solid 




Figure 10: (Color on line) Average shear stress at the inner cylin- 
der S as a function of the confining pressure P for different cylinder 
roughness: (I, O) sandpaper, (*) glued seeds and (k, A) corrugated. 
The full symbols indicate smooth horizontal walls, while hollow ones 
indicate rough horizontal walls (Vg = 3.14 mm/s) 



fraction changes. This indicates that most of the vol- 
umetric variations are localized in the region near the 
cylinder, as reported by other experiments and simula- 
tions SSlm 



5. Conclusions 

As a way to study the rheology of granular materials 
near an interface, we have developed a cylinder shear 
apparatus inserted in the MRJ device, which allows the 
control of the radial confining pressure exerted by the 
outer wall on the grains and the measurement of the 
torque on the inner shearing cylinder. We have also 
developed a multi-layer velocimetry sequence for the 
simultaneous measurement of velocity profiles in dif- 
ferent sample zones. This study shows that the rough- 
ness of the shearing interface significantly affects the 
granular interface rheology (sliding velocity and fric- 
tion coefficient), while the influence of the roughness of 
the transverse confining walls remains localized in the 
very first layers. Moreover the measurement of the solid 
fraction profiles, based on static imaging of the sample, 
shows that, when starting to shear the material, its volu- 
metric variations remains localized close to the shearing 
interface. Those data are useful to test models attempt- 
ing to describe granular interface rheology 11641 Il5l 16511 
for annular shear. Using this device together with the 
velocimetry sequence, work is under progress to study 
the transient behavior of the granular material, under 




Figure 11: Static images used for the measurement of the solid frac- 
tion profile (sample quarter). Definition: 25 pixels/mm 2 , (a) Before 
and (b) after shear (corrugated cylinder). 



monotonic or cyclic shear, and the influence of liquid 
saturation. 
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